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H We construct a number of lattice fermions, which fulfill the Ginsparg- Wilson relation 

either exactly or approximately, and test them in the framework of the 2-flavor Schwinger 
model. We start from explicit approximations within a short range, and study this for- 
mulation, as well as its correction to an exact Ginsparg- Wilson fermion by the "overlap 
formula" . Then we suggest a new method to realize this correction perturbatively, without 
using the tedious square root operator. In this way we combine many favorable properties: 
good chiral behavior, small mass renormalization, excellent scaling and rotational invari- 
ance, as well as a relatively modest computational effort, which makes such formulations 
most attractive for QCD. 
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1 Introduction 

Recently there have been intensive new activities to circumvent the notorious Nielsen- 
Ninomiya No-Go theorem ||I| for chiral fermions on the lattice. That theorem has a rather 
complicated set of minimal assumptions; we simplify them slightly to the following state- 
ment: for an (undoubled) lattice fermion, unitarity, discrete translation invariance, locality 
and (full) chiral symmetry cannot coexist. To get around this theorem, the breaking of 
each of these properties has been tried. Examples are (referring to the above order) one- 
sided lattice Dirac operators, random lattices, SLAG and Rebbi fermions, and finally the 
Wilson fermion. 

The latter breaks the full chiral symmetry in a rather hard way, destroying essential 
physical properties related to chirality. (By full chiral symmetry we mean the relation 
{D,75} = 0, where D is the lattice Dirac operator, and the curly bracket denotes the anti- 
commutator.) The subject of this paper is a recently re-discovered approach to perform 
such a chiral symmetry breaking in a much softer way, preserving a modified but continuous 
form of chiral symmetry at finite lattice spacing. 

That remnant lattice chiral symmetry transformation has first been written down in 
Ref. 0. Its generalized form reads (we use a short-hand notation for the convolutions in 
coordinate space) 

i,^^(i,{l + e[l-DR]^^))^, ^, ^ [{I + e[l - DR]-f^) ^) ^ (1.1) 

where i? is a local Dirac scalar, i.e. it decays at least exponentially. We now require 
invariance of the fermionic Lagrangian ipDip to 0(e). This amounts to the condition 

{L»,,j„75} = 2(L»75i?I?).,y , (1.2) 

which is known as the Ginsparg- Wilson relation (GWR). The crucial point is that the 
term R, which describes the chiral symmetry breaking of D^^, 

Rx,y = \l^{D-l,^5} . (1.3) 

is local, which is not the case for Wilson fermions, massive fermions etc. Therefore the 
pole structure of D~^ is not affected by R. 

Ginsparg and Wilson pointed out a long time ago that this is a particularly soft way to 
break chiral symmetry on the lattice Q . Renewed interest was attracted to this approach 
especially by Ref. 0]. It was demonstrated that it preserves the triangle anomaly |p, 
that it avoids additive mass renormalization and mixing of matrix elements [^ , and that 
it reproduces the soft pion theorems ||6|. So far we refer to vector theories, but the 
GWR even serves as a basis for the construction of chiral gauge theories [0]. As further 
applications, the GWR provides for instance a safe continuum limit of the chiral anomaly 
Ig] and of the spontaneous chiral symmetry breaking Q. It also opens new perspectives 
in other fields, like random matrix models pO|]. 



As a general ansatz for a solution of the GWR, we write 

D-^ = D-^ + R, where {L'-\ 75} = . (1.4) 

Now a suitable operator D^ — with the correct continuum limit — has to be identified, 
which is non-trivial even in the free case. If it suffers from doubling, non-unitarity or 
SLAG type non-locality (finite gaps of the free Dy^{p), where p is the momentum), then 
that disease is inherited by D. Hence all these options must be discarded. A way out is. 



however, a non-locality of the Rebbi type |11|, where D^(p) has divergences. This still 
allows for locality of D. 

That is exactly the mechanism which is at work to yield local perfect fermions [^ 
(lattice fermions without lattice artifacts). Here the term R~^ plays a specific role: in 
the block variable renormalization group transformation, which leads asymptotically to a 
perfect action, it is the kernel between the blocks in a Gaussian transformation term. R 
Also in the interacting case the perfect fermion is a solution of the GWR, in agreement 
with the fact that it breaks the chiral symmetry only superficially (in the manifest form 



of the lattice action), but not with respect to the physical observables [12|. However, the 



perfect fermion can only be constructed perturbatively or in the classical approximation 



(fixed point action, FPA) ||T^, 13, 14|. It turned out that the latter is a solution of the 



GWR too [Ql, and its locality is optimized by the choice 

Rx,y = -Sx,y , {^x,j/,75} = {D-f^D)^^y . (1.5) 

The use of a FPA, together with the corresponding classically perfect topological charge, 
guarantees that the index theorem is correctly represented on the lattice 0]. 

We refer to eq. ( p..5| ) as the ^^standard forrrH^ of the GWR. For a resulting FPA, the 
index theorem has been confirmed numerically in the Schwinger model |15], where the 
topological charge was also defined in the spirit of the FPA. 

Except for the FPA, another GW fermion was discovered by H. Neuberger [p^]. The 
locality of that solution has been established analytically in a smooth gauge background 



and numerically up to moderate coupling strength in QGD |17|. 



The relaxation of the full chiral symmetry condition to the GWR allows the fermions 



to be local, but it has been conjectured that they can still not be "ultralocal" |15]. This 
means that their couplings may decay exponentially, but they cannot stop at a finite 
number of lattice spacings, not even in the free case. In fact, this has been demonstrated 
for the standard form of the GWR jl^ and for a more general class of GW kernels R in 
any dimension d>2 [pO| ]. 

Against this background, we construct in the next Section lattice fermions, which 
satisfy the GWR at least to a good approximation, but which have couplings only in a 
short range, so that they can be implemented directly for simulations. ^ At the same 

^Only the limit i? — > ("5 function blocking") leads to non- locality. 

^ There is some conceptual — though not technical — similarity with the papers in Ref. |2l[. 



time, we optimized other essential properties, in particular the scaling behavior and the 
approximate rotation invariance. We emphasize that the GWR does not guarantee a good 
quality of those properties. In fact, it has been observed p^, |2^ — and it will be confirmed 
in Section 3 — that the Neuberger fermion does a rather poor job with that respect. 

In addition we do our best to keep the computational effort modest. In the framework 
of the Schwinger model, where our study takes place, we can also simulate very complicated 
actions. However, we are interested in working out a formulation, which has the potential 
to be carried on and applied to d = 4. All the above properties (in particular excellent 



chirality and scaling) are also manifest in the FPA for the Schwinger model in Ref. |14], 
except for the simplicity. That (approximate) FPA involves 123 independent couplings — 
and all together 429 terms per site. Unfortunately this makes an analogous formulation in 
d = 4 inapplicable. So far, all attempts to construct a useful approximate FPA for QCD 
got stuck in the fermionic part of the action, hence it is strongly motivated to search for 
a simplified alternative with similar qualities. 

2 Hypercubic approximate Ginsparg-Wilson 
fermions 

To improve the lattice fermion, i.e. to suppress lattice artifacts, we have to include cou- 
plings to lattice sites beyond nearest neighbors. However, the systematic extension to 
various lattice spacings amplifies the number of couplings very rapidly. As an option, 
which seems to allow for a powerful improvement, but which is still tractable in QCD 
simulations |24, 25], we focus on the ^^hypercube fermion!^ (HF), where Tpx is coupled to 



tpy ii \Xfj_ — y^\ < 1 for // = 1 . . . d, i.e. we couple all sites inside a d dimensional unit 
hypercube. 

2.1 Free hypercube fermions 

For free fermions, perfect actions can be computed analytically [13|. However, they can 
only be local in the sense that their couplings decay exponentially, in agreement with 
the conjecture about the absence of ultralocal GW fermions. For practical purposes the 
parameters in the renormalization group transformation can be tuned so that the decay 
becomes very fast. This is the case for the parameters corresponding to the standard 
GWR. There the couplings were truncated to the unit hypercube by means of periodic 
boundary conditions over three lattice spacings p6| . The resulting truncated perfect HF 
(TP-HF) has strongly improved scaling properties compared to the Wilson fermion. Since 
the truncation is only a small modification, it also approximates the standard GWR to a 
good accuracy ]l8| . 

To fix our notation, we write the free lattice Dirac operator as 

D,^,+r = D{r) = /5^(r)7^ + A(r) , (x, r G Z'^) , (2.1) 





TP-HF 


CO-HF 


SO-HF 


7^ 


= Pi(l,0) 


0.30938846 


0.30583220 


0.334 


p(^) 


= />i(l,l) 


0.09530577 


0.09708390 


0.083 


Ao 


= A(0,0) 


1.48954496 


1.49090692 


1.5 


Ai 


= A(1,0) 


-0.24477248 


-0.24771369 


-0.25 


A2 


= A(1,1) 


-0.12761376 


-0.12501304 


-0.125 


Vst 


3.008 • 10-4 


1.006 • 10-4 


65.518 • 10-4 



Table 1: The couplings of free hypercube ferniions — truncated perfect, chirally optimized 
and scaling optimized — and their violations of the standard GWR. Note the constraints 



p^ ' + 2p^ ' = 1/2 and Aq + 4(Ai + A2) = from the continuum limit and from mass zero, 
respectively. 



where we assume the sensible symmetries: p^ is odd in the p direction and even in all 
other directions, while the Dirac scalar A is entirely even. In addition p^ is invariant under 
permutations of the non-p axes, and A under any permutation of the axes. Of course D 
must also have the correct continuum limit, D{p) = ip^^^ + 0{p'^). 

As a measure for the total violation of the free GWR, we sum the squared violations 
in each site. 



V = J2[Y.'^D(^h5R{y)D{z)6r,^ 



+y+z 



{D{r),l^] 



r x,y,z 



4E [ E [A(x)i?(2/)A(z) - p^{x)R{y)p^iz)]5r,,- 



x,y,z 



-y+z 



A(r) 



For the standard GWR (p..5|) this simplifies to 



"^«* = E [E['^(^)^(^ -x)- p^{x)pi,{r - x)] - A(r) 



(2.2) 



(2.3) 



If we just optimize the couplings so that Vst becomes minimal, then we can still do some- 
what better than the TP-HF. For d = 2 this can be seen from Table |l[ where we denote 
that ^^ chirally optimized" HF as CO-HF. However, the chiral optimization makes the scal- 
ing behavior a little worse than it is the case for the TP-HF, see below. On the other 
hand, if we search for excellent scaling, then we end up with a HF that we call SO-HF 
{scaling optimized), and which is also included in Table ||. Its value for Vst is still small, 
but clearly larger than in the previous two cases. The TP-HF can therefore be seen as a 
compromise between the two optimizations with respect to just one property. Note that 
the couplings are similar in all these three HFs. 

Actually the comparison of Vst is not really fair for the SO-HF, because that fermion is 
not necessarily related to the standard GWR. If we allow for R^^y = r^dx^y (ro arbitrary), 



then its value V drops to 18.413 • 10"^ (for tq = 0.48950632). If R is generalized further 
to an even hypercubic form (like A), we arrive at V = 15.977 • 10~^, but this is still not in 
the same order of magnitude as the TP-HF and CO-HF. 

We can also perform the minimization of the HF couplings for such generalized GW 
kernels R. However, we have to avoid the trivial minimum at A = i? = (naive fermion), 
so we insert another constraint on A. We require the mapping to d = 1 to reproduce 
the Id Wilson fermion, Ai + 2A2 = —1/2 (which is also the case for TP-HF). The min- 
imum is still found in the same vicinity, p^"^' = 0.09909187, A2 = —0.13035191, and 
/?(0,0) = 0.50514453, i2(l,0) = -0.00141375, R{l,l) = 0.00183004, which shows that — 
at least with this extra constraint — the vicinity of the standard GWR plays indeed a 
special role. The violation then amounts to V = 0.793 • 10"'*, so we can apparently not 
proceed to lower orders of magnitude any more. 

Hence we return to the form Rx,y = r^dx^y In this case, there is a simple way to 
illustrate the accuracy of the GWR. For an exact GW fermion, the spectrum lies on 
the circle in the complex plane with center and radius l/(2ro) [p. This holds with or 
without gauge interaction. Therefore we can just check how close the eigenvalues of the 
approximate GW fermion are to that circle. This is shown in Fig. ||, which confirms the 
hierarchy measured from V. The strongest deviations occur in the arc opposite to zero. 
This arc corresponds to high momenta and therefore a fine resolution, which is sensitive 
to small inaccuracies of the GWR. 




TP-HF 




CO-HF 



SO-HF 



Figure 1: The free spectra of 3 hypercube fermions: truncated perfect (left), chirally 
optimized (center) and scaling optimized (right), on a 100 x 100 lattice. 



However, as we emphasized in the introduction, we do not want to concentrate solely 
on the chiral quality. In particular we want to consider the scaling behavior as well, and 
we now take a first look at the impact of such a GWR optimization with that respect. 
Fig. H shows the dispersion relations for the free HFs mentioned before, and we see that 
all the three are strongly improved compared to the Wilson fermion. We also see that the 
hierarchy among the HFs is inverted. This observation motivated the consideration of the 
SO-HF, which was constructed by hand. (Some comments on the construction are given 
in the appendix.) 



continuum 

Wilson 

TP-HF 

CO-HF 

SO-HF 



LU 




3.5 

3 
2.5 

2 
1.5 

1 

0.5 




Figure 2: The dispersion relation for free hypercube fermions compared to the Wilson 
ferniion and to the continuum. (The energies of all higher branches keep above 4.4.) 



We also want to compare thermodynamic scaling properties of the free HFs. In Fig. 
y we show the ratios P/T"^ at /Xc = 0, resp. P/ I-l^ at T = 0, which are scaling quantities 
in d = 2 [P: pressure, T: temperature, ^c- chemical potential). In the HF actions, fic is 
incorporated according to the prescription in Ref. [^3- For T — > (many lattice points 
Nt in Euclidean time direction), resp. ^c — > 0, all HFs converge to the continuum ratio, 
but the speed of convergence differs strongly. The relative quality of the scaling behavior, 
which was suggested by the dispersion, is confirmed. Again the SO-HF looks particularly 
impressive. 

2.2 Applications to the Schwinger model 

We now proceed to the 2-flavor Schwinger model (QED2), and we attach the free fermion 
couplings equally to the shortest lattice paths only. Moreover we add a clover term, which 
turned out to be useful, and we fix its coefficient to csw = 1- 

The classically perfect fermion-gauge vertex function includes a number of plaquette 
couplings spread over some range |2^. However, if we sum them up and "compactify" 
them all into the clover term, we obtain csw = 1- This holds for the "compactified" 
on-shell 0{a) improvement of any of our massless HFs. Note that — unlike QCD — this 
value does not get renormalized due to the super-renormalizability of the Schwinger model 
|28|. Therefore, csw = 1 provides a non-perturbative 0{a) improvement for the Wilson 
fermion as well as the HFs. 

We use quenched configurations on a 16 x 16 lattice. However, in the evaluation of 
dispersion relations and correlation functions (see below) the square of the determinant is 
included as a weight factor, following the prescription in Ref. |22|, Section 3. For the pure 
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Figure 3: Thermodynamic scaling ratios. Left: The ratio P/T'^ (at Hc = 0) as a function 
of the number Nt of lattice points in the Euclidean time direction. Right: The ratio -P//ic 
at T = as a function of the chemical potential Hc- 



gauge part, we use here and throughout this paper the Wilson plaquette action, which is 
actually perfect in 2d Abelian gauge theory [p^ . 

We first compare the spectra of our three HFs at /? = 6. Figs. ^ and |l^ show these 
spectra for a typical configuration. They are still fairly close to the unit circle (SO-HF is 
closer to a slightly larger circle, in agreement with the observation that it has an optimal 
ro < 0.5, see previous subsection). The splitting of the HF spectra around 2 can be viewed 
as a "residue" of the Wilson double circle. Close to zero we observe that the smallest real 
eigenvalue becomes finite. In QCD at weak coupling, such an effect corresponds roughly 
to the quark mass renormalization. Referring to this analogy, we denote the smallest real 
eigenvalues in the following as "quark" mass renormalization Aniq. For the HFs at /? = 6 
it amounts to ArUq ~ 0.03. We also see from Fig. |^ that the clover term has actually 
a negative impact on the eigenvalues in the region around 2, but it improves the more 
important eigenvalues close to 0. (In this context, the Sheikholeslami-Wohlert action was 
studied in the Schwinger model in Ref. |^8|; for a systematic study in d = 4, see Ref. 
||29|.) In particular it decreases Am,q, both, for the HFs and for the Wilson fermion. R 
Interestingly the mass renormalization is practically the same for the HFs and for the 
Wilson fermion with clover, and it is again almost the same if we omit the clover term, as 
it was observed before in QCD [pqj . 

Of course, the mass renormalization — and the deviation from the circle in general 
— increases at stronger coupling; as an example we show a typical TP-HF spectrum at 
/3 = 4 and at /3 = 2 in Fig. |. 

As a scaling test, we consider the dispersion relations of the "mesons"; we call them 
"vr" (massless) and "r/" (massive). In Fig. ^ we show the vr and r/ dispersions for our 

^Here we refer to the non-critical Wilson fermion: p'^' — — Ai — 0.5, p^'^' = A2 = 0. 
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Figure 4: Spectra of different HFs for a typical configuration at /? = 6, compared to the 
Wilson fermion with and without clover term. (For the SO-HF with the same configu- 
ration, see Fig. |I^. j The clover term is always included in the HF actions if not stated 
differently. It has a negative impact in the region around 2, but it improves the more 
important region close to 0. 



three HFs, and we compare them to the results for the Wilson fermion and for the FPA. 
For the Wilson fermion we now use the critical hopping parameter Kc = 0.25927, which 
was determined using the PCAC relation [^]. These dispersions were obtained from 5000 
configurations, using the same ensemble for all the fermion types in Fig. |^ (as well as Fig. 
11 in Section 3). We see that our HFs are all strongly improved, to the same level as the 
FPA. This is very remarkable, because we only use 6 independent terms per site — as 
opposed to 123 in the FPA. (By contrast, the dispersions for the Sheikholeslami-Wohlert 
action are very similar to the Wilson action, see first Ref. in |2^.) In addition to the 
excellent scaling — in particular for the SO-HF — the r] dispersion also reveals a good 
agreement with asymptotic scaling, which predicts an ij mass of ?n,^ = ^2/(7r/3) ~ 0.326. 
As a further test for the quality of our fermion actions, we show in Fig. ^ the decay of 




0.5 1 1.5 



Figure 5: Typical spectra of the TP-HF at P = 4 and (3 = 2. The resuhs for the CO-HF 
and SO-HF are similar. Strong coupling amplifies the deviation form the unit circle, and 
in particular the mass renormalization. 



the correlation function 



C3(x) = (Vi(0)a3^(0)-V^(x)a3V(2;)). 



(2.4) 



We are particularly interested how well rotational invariance is approximated (which can 
only be measured by using o"3 |jl^ ) . All our HFs yield a very smooth decay — again on the 
same level as the FPA — which reveals an excellent approximation of rotational invariance 
already at short distances. (The jumps around |x| ~ 8 are finite size effects). On the 
other hand, for the Wilson fermion the decay performs a zigzag at short distances, which 
corresponds to a "taxi driver metrics" q it takes rather large distances to approximate 
the Euclidean metrics instead. 

Finally we also compared the percentage of configurations, which fulfill the index the- 
orem, if we use the geometric definition for the topological charge. Here the results for the 
SO-HF are similar to those for the Wilson fermion at /? = 6, 4 and 2, which were reported 
in Ref. P]. The percentages amount to 100 %, 99.5(2) % and 78.7(13) %, respectively. 



Our HFs are successful with respect to scaling and rotational invariance, and they 
approximate the GWR reasonably well up to moderate coupling strength. Their one 
unpleasant feature is the quite significant mass renormalization, which is visible again in 
Fig. |6|. For the n mass it amounts to m-,^ ~ 0.13 (which is consistent with Anig). This is 
a practical problem: in this formulation, ni^^ = would require the tuning of a negative 
bare mass. For the massive TP-HF that is also unfavorable for locality. However, we do 
know the TP-HF at any finite bare mass, and we can construct also massive SO-HFs, see 

^In the "taxi driver metrics", the distance between two lattice sites is given by the shortest lattice 
path(s) connecting them. 
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Figure 6: Meson dispersion relations at fi = 6: Wilson fermion (diamonds), TP-HF (empty 
circles), CO-HF (triangles) and SO-HF (little boxes) — all the HFs with a clover term — 
compared to the FPA (filled circles) and the continuum (solid line). 

appendix. Here we do not simulate them — that would also contradict our attempt to 
stay close to the GWR — but in d = 4 this option should be reconsidered in order to 
tackle the problems related to the mass renormalization, which are more difficult there. 

Indeed, the same problem was a major obstacle in similar QCD simulations \2G, 24, 31, 
25|. There the critical bare mass could be shifted to by the use of fat links with negative 



staple terms |32|. We also performed tests with fat links in the present framework, using 
both, positive and negative staple terms, and we could remove the mass renormalization 
here too. It is profitable to distinguish a staple term with weight wiXi attached to the 
scalar term, and an independent staple weight WrP^^' attached to the vector term. Then 
the weight of the direct link is (1 — 2wi)Xi resp. (1 — 2wi)p^^> . In Table g we give 
results from various staple combinations for the "quark" mass renormalization Atti^ and 

® An attempt to optimize the fat links analytically — so that the GWR violation is minimized — yielded 
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Figure 7: Tije decay of the correlation function |C3(x)| (defined in eq. \2.^) , illustrating 
the level of approximate rotation invariance for various fermion actions. 



for the mean value of the squared distance from the unit circle, {5f.). The first quantity 
refers to the (physically essential) small real eigenvalues, and the second quantity to the 
whole spectrum. In particular, Aniq = requires a strongly negative wi, which amplifies, 
however, ((5^). As a remedy for that we can choose a positive Wr (and modify wi a little). 
However, here we optimize the GWR only, and it turned out that "GWR optimal" staples 
are unfavorable for the scaling behavior. The vr dispersion relation looks still fine, but the 
r] dispersion is somewhat distorted. For that reason — and for the sake of the simplicity of 
the action — we do not use fat links in the rest of this paper, i.e. we return to wi = Wr = 0. 
Instead, we are going to suggest other solutions for the mass renormalization problem 
in the following two Sections. The goal is to further reduce the GWR violation, since Arrig 
and rriTT vanish continuously as V approaches 0. This is in agreement with the consideration 
in Ref. Q, where a small fermion mass is introduced as a regularization before taking the 
chiral limit (in the sense of the GWR) . 

3 Exactly massless fermions with a good scaling behavior 

Let us introduce the operator 



V = l--D, ifi^O) 



(3.1) 



positive wi, Wr, which led, however, to even larger values for Atti, 
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Wl 


Wr 


Arriq 


{6'^) (in units of 10"^) 








0.032(8) 


1.34(4) 





0.1 


0.029(8) 


1.07(3) 





0.2 


0.027(7) 


0.91(3) 


-0.1 





0.021(6) 


1.46(5) 


-0.3 





-0.001(4) 


2.07(14) 


-0.24 


0.2 


0.002(4) 


1.59(11) 



Table 2: Two characteristic quantities to measure the approximation of the standard 
GWR: the "quark" mass renormahzation Am,q, and the mean value of the squared radial 
deviations of the eigenvalues from the unit circle (in the complex plane), ((5^). We show 
results for varying fat links, with the scalar staple coefRcient wi and the vector staple 
coefficient Wr, at /3 = 6. Criticality can be achieved by a strongly negative wi. 



where /x is a real mass parameter, and we require the GWR to hold with Rx^y 
This is equivalent to 

The overlap type of solution inserts an operator V of the form 



V 



A 



\fWA 



(3.2) 



(3.3) 



If this is well-defined, then V is unitary and the GWR reads V^ 
lattice Dirac operator Dq as 

A = ^x-Do , 



75^75- Inserting some 
(3.4) 



leads to a correctly normalized Dirac operator D. Most of the literature deals with the 
Neuberger fermion, which uses Dq = D\y, where Dy^/ is the Wilson-Dirac operator. Neu- 
berger fermions were simulated (quenched) in the Schwinger model p^ , p^ , 23] and in 4d 
non-Abelian gauge theory ||3^, |l^. The mass parameter // is usually set to 1. As long as 
we stay with the Wilson-Dirac operator it can be absorbed in the mass M of D^. But if 
we insert a more sophisticated Dq, then fi takes a non-trivial role, see below. 

It was observed that the use of a suitable Dq = Dhf instead of Dw improves the 
locality of the free D significantly [18|, and it looks promising also with other respects. 
Further properties of improved overlap fermions will be discussed in this Section. 
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3.1 Free overlap fermions 

If we insert the ansatz Dq = p^j^ + A, then the free overlap Dirac operator and its inverse 
read (in momentum space) 

^P^ + V-/^^(p) + [A(p)-^P ' 

D-Hp) = K^"Q§H^)' Q(P) ■■= \/-p'(P) + ^'(P)' ^(P) ■■= ^(P) - /^ • (3-5) 

There are two conditions for a pole in the free propagator: 

(1) p\p) = 0, (2) u{p)<0. (3.6) 

Hence the vector term p^ determines the shape of the fermion dispersion, and the scalar 
term A can fix an end-point of that dispersion curve [|. (This is different from the usual 
case, where both terms contribute to the shape, and the curve cannot just end inside 
the Brillouin zone.) From the symmetry properties that we assumed for the vector term, 
condition (1) allows for poles whenever all momentum components obey p^ G {0, tt}, hence 
it allows for fermion doubling with 2 species. However, the condition (2) can still save 
us from doubling; A and fi should be chosen such that only the pole at p = really 
occurs. If we insert D\y with a mass M, then the condition for one free species is simply 
< {fi — M) < 2. In gauge theory, some tuning of /j, might be required to stay with one 
species Q. 

For a HF with an even scalar term A, the conditions read fj 

P-Xo- 4(Ai + As) > , p - Ao + 4A2 < , p - Aq - 8A2 < , (3.7) 

where we use the notation introduced in Table ||. In the presence of gauge interactions, 
the 1-species region becomes more narrow from both sides. This can be understood from 
a crude consideration, which introduces "effective links" somewhat smaller than 1, and 
which therefore moves the (negative) couplings Ai ~ 2A2 to larger "effective" values. This 
was also observed numerically for the TP-HF in the Schwinger model ]3^. q Hence we 
should start far from the boundaries in the free case, to be on more solid grounds in gauge 
theory. The truncation by condition (2) should occur not too close to the edge of the 
Brillouin zone, but of course quite far from p = 0, because this region is really needed. 
It appears that // = 1 is a very reasonable choice. In Fig. ^ we show the free fermion 
dispersion for overlap fermion D constructed from Dw and from our three variants of 
Dhf (all of them massless): we show the full curve given by condition (1), and we mark 
the end-points for p = 1. Those end-points can be shifted arbitrarily by changing /x. 
(Higher branches do not appear, see appendix.) 

'^It is the square root which causes this unusual behavior. 
*The situation on the boundaries is tricky. 

'At this point, we do not insist on ^ A(r) — 0, i.e. Dhf rnay also be massive. 

^"As an example, the massless TP-HF has the free condition is < /i < 1.96. At /3 = 3 the physical 
interval shrinks to about /i £ (0.3, 1.8); outside this interval the index theorem is violated. 
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Figure 8: The free dispersion relation for the Neuberger fermion (inserting Dq = D]y in 
the overlap formula) and for various improved overlap fermion (inserting Dq = Dhf)- We 
mark the end-points of the curves for the mass parameter // = 1. 

We see that the overlap dispersion based on Dw is unfortunately worse than the 
dispersion of the ordinary Wilson fermion, cf. Fig. |^. (Technically, the reason is that in 
the ordinary case the scalar term helps to raise the curve a little). It has been conjectured 
that overlap ferniions might be free of 0{a) artifacts even in gauge theory. However, for 
the Neuberger fermion the 0{a'^) artifacts are quite bad, as we see even in the free case. 
The improvement by using Dhf is again very significant, in particular for the SO-HF. 

The improved scaling of free overlap-HFs is also confirmed if we repeat our thermody- 
namic considerations, see Fig. p. The good quality of these fermions does hardly come as 
a surprise [^] : if we insert a fermion into the overlap formula, which obeys the GWR for 
R = '^a-,j//(2/i), then it reproduces itself due to A^A = fi^. Now we insert an approximate 
GW fermion, so its modification by the overlap formula is rather modest, and the good 
scaling quality does essentially persist, [j 

In the appendix we comment on the use of massive Dhf — which still produce massless 
overlap fermions — and also on the option to render overlap fermions massive themselves. 

3.2 Interacting overlap ferniions 

Also in the interacting case, the overlap-Dirac operator 

A 



D 



fJ- 



1 



VA^ 



(3.8) 



"^Here the overlap correction A -^ fiA/y'AKA is completely obvious. 
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Figure 9: Thermodynamic scaling ratios: P/T"^ at fie = (left), and P/fJ,'^ at T = 
(right), as in Fig. y, but now for the corresponding overlap fermions. Again the overlap 
SO-HF is clearly most successful. 



satisfies the GWR with Rx,y = '5x,y/(2/i), hence its spectrum is situated on the circle with 
center and radius l//i. We now wonder what really happens to the eigenvalues of Dhf 



as the hypercube fermion is inserted in the overlap formula. Fig. 10 shows as an example 
those eigenvalues (for a typical configuration at /3 = 6 on our 16 x 16 lattice) before and 
after application of the overlap formula. Especially in the region around zero — where the 
eigenvalue density is low — there is an obvious mapping of the eigenvalues one by one onto 
the circle. In such cases, the effect of the overlap formula is close to a radial projection of 
each eigenvalue onto the circle. This can easily be understood from the overlap formula, 
and it is in agreement with the above statement that the overlap modifies an approximate 
GW fermion only modestly. (Of course this observation is not so obvious if we insert a 
fermion far from a GW fermion, such as D^y.) 

This geometric understanding of the overlap formula also provides a neat interpretation 
of the parameter fi: it is simply the center of the circle (through 0) that the spectrum is 
projected on. Of course that center must be chosen between the small real eigenvalues of 
the original spectrum (which have to be mapped to 0) and the larger real eigenvalues (to 
be mapped on 2fi). So the allowed range of // can be recognized immediately from the 
spectrum of Dq, see e.g. Figs. ||, |^, ^ and IC. This range shrinks at stronger coupling. If 
we do not respect it, then we map one or several eigenvalues to the wrong arc crossing the 
real axis, and as a consequence the index theorem is violated (as we mentioned already in 
the previous subsection). 

These observations raise hope that also in gauge theory the scaling of an approximate 



GW fermion is not much affected by the ^^chiral proiectiori\ Indeed, Figs. |ll| and |l2|show 
that the scaling quality — tested by the meson dispersion — as well as the approximate 
rotational invariance is still very good, whereas the Neuberger fermion is contaminated 
by considerable artifacts (again it is a little worse than the ordinary Wilson fermion). In 
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Figure 10: The mapping of the eigenvalues of the SO-HF (for a configuration at f3 = 6) 
onto the unit circle (in the complex plane) by means of the overlap formula. 

particular for the overlap SO-HF the vr dispersion is still excellent, but the tj dispersion 
is a little distorted by the overlap formula (the slope is a bit too steep). The same effect 
occurs if we remove the mass renormalization by means fat links (using a negative wi , cf . 
Section 2.2). 



As a further issue, we now want to address the question of locality. For the Neuberger 
fermion in d = 4, locality was established analytically for very smooth gauge fields |17]. 
At fi — M = 1 the dimensionally generalized condition is that any plaquette variable P 
has to obey 



\l-P\\ < 



(3.9) 



5d{d - 1) ' 

since this rules out the danger of A'^A = 0. We see that this condition is more powerful 
in d = 2 than in d = 4. 

In numerical tests in QCD down to /3 = 6, it turned out that the degree of locality is 
only reduced gradually by the gauge interaction [l^. R The observation that it remains 
rather close to the free locality raises hope that the improvement of free locality persists 
also in gauge theory. In fact, this is confirmed in our Schwinger model study as we see 
from Fig. |l^. It compares first the free locality, and then the absolute value of the maximal 
correlation over a distance \x\ at /3 = 6 on a 24 x 24 lattice. More precisely, we show the 



As the interaction is turned on, the value of fi 
above 1. 



M, which is optimal for locality, moves somewhat 
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Figure 11: Meson dispersion relations at (3 = 6; Neuherger fermion (diamonds), overlap 
TP-HF (open circles) and the overlap SO-HF (little boxes) compared to the continuum 
(solid line). 



expectation value of 

f(r)= Ti 11^^(^)11 I l^-yl=r} (3.10) 

for a unit source at x, as it was done before for the Neuberger fermion in QCD [|l^. We 
see that the exponential decay is much faster for the overlap SO-HF than it is the case for 
the Neuberger fermion. (Again, the httle bumps in the middle of that plot are finite size 
effects; they move continuously if we vary the lattice size.) 

A possible danger for the locality of an overlap operator could still occur at strong 
coupling, if the eigenvalues of A'^A cluster very densely close to zero. Since the use of an 
exact GW fermion (with respect to Rx,y = Sx,y/{2fi)) for Dq fixes A^A = /i^ = const, (for 
any configuration), one could have hoped that an approximate GW fermion suppresses the 
density of eigenvalues near zero. However, from the eigenvalue histograms we could not 
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Figure 12: The decay of the correlation function |C3(x)| (defined in eq. 2A), iUustrating 
the level of approximate rotation invariance for two overlap fermions. 



find any significant difference between tlie Neuberger fermion and tlie overlap SO-HF with 
this respect at /3 = 6, and at stronger coupling neither |35]. In Fig. ^ we show histograms 
of eigenvalues obtained from 1000 configurations at /3 = 6: the total eigenvalue density 
(left) and the lowest two eigenvalues (right). Also if we measure the separation of those 
lowest two eigenvalues in each configuration, there is no significant difference between the 
Neuberger fermion and the overlap SO-HF. (Apparently the different degree of locality 
can be recognized from the small eigenvalue distribution only if one takes into account the 
(topological) quality of the considered eigenvalues.) 

As a last comparison, we show the angular density of the overlap SO-HP and the Neu- 
berger fermion (at /3 = 6, using fj, = 1) compared to the FPA in Fig. |l5|. qj Since there 
is a refiection symmetry on the real axis, we only show the upper part, ranging from 
{6 = 0) to 2 {9 = it). At small angles, i.e. small momenta, there is hardly any difference. 
Again a clear distinction occurs at the opposite arc. The overlap SO-HF has a peak at 
practically the same angle as the fixed point fermion, in contrast to the Neuberger fermion. 

To summarize this Section, we first repeat that the overlap-HFs are exact GW fermions. 
If they scale very well in the straight application and they approximate the GWR decently 
— like the SO-HF — then the resulting overlap fermions do still scale well and they have 
an exact remnant chiral symmetry. So they combine two very important properties, which 
were simultaneously present only in the FPA so far. It's draw-back, however, is that its 
simulation is still quite involved — although the action is much simpler than the FPA. 
In d = 4 the square root operator is generally problematic, and in addition we use a 



For the Neuberger fermion, such an angular density was studied before in Ref. 
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Figure 13: The locality of the overlap SO-HF compared to the Neuberger fermion. Left: 



the decay of the free couplings (as introduced in eq. (2A)). Right: the decay of the 
maximal correlation (as defined in eq. ( ]3.i0|) ) at f3 = 6 on a 24 x 24 lattice. 



matrix A, which is not as sparse as it is the case for the Neuberger fermion. Therefore the 
apphcation of such a formulation in QCD might be rather expensive. In the next Section 
we present a suggestion on how to reduce the computational effort. 

4 Perturbative chiral correction 



We have seen in Section 2 that suitable HFs can scale well and approximate the GWR 
reasonably well, up to a certain coupling strength. In Section 3 we converted these HFs 
into exact GW fermions by means of the overlap formula, without a strong distortion of the 
good scaling and rotation invariance. However, this formulation contains an inconvenient 
square root operator. We now suggest a new method to avoid that operator in order to 
reduce the computational effort. For that purpose, we perform the "chiral projection" 
perturbatively. 

Assume that we start from an operator Dq, which is close to a GW fermion for Rx^y = 
6x,y/{2fi). Then the operator 

e:=A^A-fi'^ (4.1) 

is small, ||e|| <C 1, and we use it as an expansion term to approximate {A^A)'^'"^. The 
perturbative chiral projection takes the form 



0{e) 



H-AY 
1 
2 



Y 



:A^A 



/"^ 



15 



'^A^A^^^AUf 



etc. 



(4.2) 



20 



100000 



10000 



<D 
CO 

> 

<D 
D) 

O 




1000 r 



100 



1000 



100 



-- 10 



0.5 



Figure 14: The eigenvalue distribution of A' A at (3 = Q for the Neuberger fermion (bold) 
and the overlap SO-HF (thin) from 1000 configurations. Left: all eigenvalues. Right: the 
lowest two eigenvalues (solid and dashed). 
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Figure 15: The angular eigenvalue distribution on the unit circle for various GW fermions, 
measured in 5000 quenched configurations. 

One could write down the fully explicit action corresponding to the Dirac operator Dpcp, 
but this form is complicated already for 0{e) (at least the absence of fat links is profitable 
here). However, for practical applications such an explicit form is not needed, hence we do 
not write it down here. Once the matrix-vector products Ax and A'x are implemented, 
the implementation of the n^^ order chiral correction is trivial, and requires essentially 
1 + 2n such multiplications. Taking the first few orders is still much cheaper than the 
exact chiral projection discussed in Section 3. We emphasize that the linear growth in n 
is very modest, so it should be feasible to simulate this expansion also beyond the lowest 
orders. 

Let us now discuss the efficiency of the perturbative chiral projection. We consider the 
SO-HF, which scales very well, but which violates the GWR most (among our HFs). As 
a first example, we consider the spectrum of the free SO-HF originally and after the first 
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order chiral projection for fj, = 1. Fig. |T^ (left) shows that this first order does most of 
the projection aheady; the resulting spectrum can hardly be distinguished from the unit 
circle. 
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Figure 16: The initial and the first order chirally corrected free spectrum (using fi = 1) 
for the SO-HF (left) and for the Wilson fermion (right) on a 30 x 30 lattice. 

Of course this method only works if we really start from an approximate GW fermion. 
P^ As an illustration, we show the spectrum of the free Wilson fermion before and after 
first order "chiral projection" in Fig. |l^ (right). Many eigenvalues do cluster at the unit 
circle (those corresponding to small eigenvalues of e), but others diverge in this inadequate 
expansion (those corresponding to eigenvalues of e with absolute values > 1). However, 
for even orders the "tail of the comet" flips far to the right-hand side, and perhaps it does 
not disturb for practical purposes. 

We proceed to the Schwinger model, and we observe that — for our HFs at /3 = 6 — 
the first order chiral correction does almost the full projection already. As an example 
we consider the SO-HF. Its initial and fully projected spectrum was shown for a typical 
configuration at /3 = 6 in Fig. |l^. We now show the first order mapping of the same 
configuration to the unit circle, and to the closest circle (of radius 1.02145) in Fig. ^ (left). 
The radius of the larger circle corresponds to the optimal value of tq , which was identified 
in Section 2.1. We observe in particular that the additive "quark" mass renormalization is 
pressed down form Aniq ~ 0.032 to about 0.014 (and to 0.002 for the second order), which 
corresponds to a reduction of rriT^ from 0.13 to 0.07 (resp. 0.02). That effect is specifically 



illustrated in Fig. 17 (right), which shows the distribution of the small real eigenvalues of 
5000 configurations initially and after the perturbative chiral correction to the first and 
to the second order. 

From Fig. |l8| (left) we see that at /3 = 4 the first order correction is still efficient, and 

^* Although simple, the expansion of the square root has not been used before, since all previous overlap 
simulations used the Neuberger fermion. 
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Figure 17: Left: The spectrum of the first order chirally corrected SO-HF for a conEgu- 
ration at f3 = 6 using fJ- = 1 (smaller circle) resp. /i = 1.02145 (larger circle). Right: A 
histogram for the small, real eigenvalues of the SO-HF in 5000 quenched configurations 
at P = 6; initial (bold), first order chiral correction using fj, = 1.02145 (solid) and the 
corresponding second order correction (dashed). The corrections move the peak — and 
hence the mass renormalization Aniq — drastically towards zero (see also Table |^. 



finally Fig. 18 (right) shows that the second order can handle even strong coupling (/3 = 2) 
quite successfully. To quantify these observations we present again the characteristic 
parameters Anig and ((5^), which were introduced in Section 2.2, see Table y. 

We have seen in Section 3 that the scaling and the approximate rotation invariance are 
not badly affected by the full chiral projection. In further tests we made the very plausible 
observation that the effect of a partial chiral projection is in between. It corresponds 
roughly to the same interpolation as we just observed in the spectrum. Since the chiral 
behavior can be controlled by this economic method, we would recommend its use most 
of all for the SO-HF, which provides very good scaling from the beginning. 



5 Conclusions 

Over the last year, the Ginsparg- Wilson relation became fashionable in the lattice com- 
munity; about 40 papers have been written about it. However, most of the literature 
focuses on the chiral properties (and recently on algorithmic questions for the Neuberger 
fermion |37, 34 1) only. Here we discussed chirality together with other crucial properties of 
fermionic lattice actions, in particular the quality of the scaling behavior, the approximate 
rotation invariance, the locality and — last but not least — computational simplicity. 
We presented three different approaches: 

• In Section 2 we discussed the straight construction of approximate GW fermions 
inside a short range. A few couplings allow for a decent approximation. Such 
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Figure 18: The spectra the SO-HF initially and after perturbative chiral correction (using 
// = 1.02145J: /3 = 4 and first order correction (left); (3 = 2 and second order correction 
(right). 

fermions are numerically tractable in QCD, and in the Schwinger model we observe 
good scaling behavior. However, at increasing coupling the GWR violation becomes 
worse, which is manifest e.g. in a stronger mass renormalization. 

• In Section 3 we corrected the GWR in the approximations of Section 2 by means of 
the overlap formula. We still observe good scaling and approximate rotation invari- 
ance for the improved overlap fermions, in contrast to the usual Neuberger fermion. 
However, this formulation is already somewhat demanding, and it is expensive to 
apply it in d = 4. 

• In Section 4 we show a way to simplify the evaluation of the overlap action, by 
avoiding the tedious square root operator. If we start from an approximate GW 
fermion, then the square root can be replaced by a simple perturbative expansion, 
and the first one or two orders are computationally relatively cheap. They do, how- 
ever, provide most of the chiral projection, i.e. the GWR violation becomes small up 
to a considerable coupling strength. At the same time, scaling and rotation invari- 
ance remain strongly improved, in particular for the "scaling optimized hypercube 
fermion" SO-HF, which appears therefore as most satisfactory. 

Based on these results, we think that an extension of this study to d = 4 is highly 
motivated. Regarding the free fermion, all our results are already relevant for d = 4 too, 
because they correspond to the special case Ps = P4 = 0. An issue in 4d gauge theory is 
for instance the choice of a suitable action for the pure gauge part. 
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Table 3: Two characteristic quantities to measure the approximation of the standard 
GWR: the "quark" mass renormahzation Am.q, and the mean value of the squared radial 
deviations of the eigenvalues from the unit circle (in the complex plane), [5'^) (as in Table 
§j. We show results for the SO-HF with chiral corrections using // = 1.02145. 
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A Improved overlap fermions with and without mass 

We start this appendix by adding some more details about the (massless) overlap fermions 
constructed from a hypercube fermion Dhf = PfiJ^i + A. If we replace the momentum 
component p2 by —lE, then condition (1) in (p.6|) yields the free dispersion relation 



coshE' 
fiPi) 



2pWp(^)fipi) ±Jl + [p(i)2 - 4p(2)2]/(p,) 



smpi 



l-4p(2)2/(pi) 
2 



with 



p(^) + 2p(^) cospi 



(A.l) 
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In all the cases considered in Section 3.1 the lower sign does not lead to a real energy E, 
hence there are no upper branches. |^ According to condition (2) this curve stops as soon 
as 

Ao + 2Ai(cospi + cosh£^) + 4A2 cospi coshE' > fj. . (A. 2) 



For Dw with mass M, this simplifies to coshi? = Wl + sin pi , ending at cospi + 

1 + sin^pi > 2 + M — yU. We observed in Section 3 that the choice /i ~ 1 sets the 
end-point to a useful position (for M = 0). 

Taking into account the normalization, the shape of the curve depends on one single 
parameter, say p^'^'. In the center of the Brillouin zone, the curve rises monotonously with 
p(^), and a good interpolation between Dw (p*-^-* = 0) and Dtp~hf (p^^-* = 0.0953) makes 
it practically coincide with the continuum dispersion up to pi ~ tt/2. In this way, the 
SO-HF was constructed, with a scalar term close to our other two HFs. We saw that it 
scales excellently, also in the straight application (without overlap). 
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Figure 19: heft: The free dispersion relation for overlap fermions constructed from massive 
TP-HFs (with mass parameter M). Again we mark the end-points for /x = 1. Right: The 
allowed region for fj, to provide one flavor for overlap TP-HFs, depending on M. 

Nothing prevents us from also inserting massive HFs in the overlap formula; we still 
obtain massless GW fermions. As an example, we use the massive TP-HF |2^, where 
the parameter M is the mass of the continuum fermion, which was originally blocked to 
the lattice. ^ If we vary M, the dispersion curve of the overlap fermion can again be 
deformed monotonously, and a value close to M = 1 appears optimal, see Fig. |l^ (left). 
The overlap fermion constructed from the TP-HF at M ~ 1 performs also very well in 
the thermodynamic tests described in Sections 2 and 3, on the same level as the SO-HF. 
However, in this region of M the 1-species range of the mass parameter // — introduced 

^^Also in all further free dispersions in this appendix, the unphysical sign never yields a real energy. 
^''For the TP-HF this implies a bare mass of M^/{exp{M) - 1). 
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in Section 3 — is very narrow, as we see from Fig. |T9| (right). |^ Since this is dangerous 
for the interacting case, we consider the SO-HF as a better approach. 

Finally we add some remarks on the case, where an overlap fermion is made massive 
by adding a mass term at the end. Such massive fermions are likely to have still very good 
chiral properties, which is useful for simulations of heavy quarks. 

We add the mass term in the straightforward manner, |^ 

Dm = l+m+ -^^ . (A.3) 

If we now search for poles in the free D^^{p), then condition (1) is generalized compared 
to (0), 

(1) p^=m{X-fi)\ rh:=l- ^ (A.4) 

{^ + m + j^Y 

while condition (2) keeps the same form. In Fig. |2^ (left) we show the free dispersion 
relations for various massive overlap fermions. (We now omit the CO-HF; its scaling is 
always slightly worse than the TP-HF.) They are all constructed by using a massless Dq 
and adding m = 1. We see that the overlap SO-HF is still good here, but not optimal 
any more. If we increase p'^' a little to 0.088, then we obtain the SOM-HF, which scales 
excellently in this case. |^ There is no reason for the optimal value of p^"^' not to depend 
on m. 

The low momentum expansion of the massive overlap operator starts with 

1 
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mo = arcosh 
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This confirms for instance the value mQ{m = 1) ~ 0.693 for all the massless fermions 
considered in Fig. 2C. In Fig. ^ (left) we plot the static vs. kinetic mass for various 



overlap fermions, and we see again that the overlap SO-HF performs well at small mo. If 



^"^The asymmetry with respect to the sign of M arises from the kernel in the Gaussian blocking term. We 
obtain optimal locality for R^^y{AI) = ±5x,yM^ /{exp{M)—M — l). Above we have chosen the positive sign. 
Then the couplings of the TP-HF at M = 1 are: p*^' = 0.11163921, p<^' = 0.02885462, Ao = 1.10520159, 
Ai = -0.09226400, A2 = -0.03854222. For the negative sign in R'^: Xi -^ -Ai. 

^* T.-W. Chin suggested an alternative way [M . In the notation of Sect ion 1 it amounts to D'^ = 



{D^+rh)/{l+RD^) . However, if _Ri:,j; oc S^.v then this is equivalent to eg. (A.3), due to Dm ~ {l+Rm)D'^, 
if m = m/{l + Rm), as S. Chandrasekharan first noticed | p5| . 

^^For completeness we give the full set of couplings of the SOM-HF: p'^' = 0.324, p'^' = 0.088; Ao = 1.5, 
Ai = -0.25, A2 = -0.125. 
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Figure 20: heft: The free dispersion relation for massive overlap fermions constructed 
from D]v or from massless HFs. Right: The free dispersion relation for massive overlap 
fermions constructed from massive HFs. In both plots we mark again the end-points for 
fj. = 1, as a useful example. 

m is of order 1, however, the two masses agree better for the overlap SOM-HF. In a next 
step, we fix the A term of the SO-HF (Ai = 2A2 = —1/4), and we tune the remaining free 
parameter p^"^' to an "ideal value", so that m^^in = "i-o- The resulting PiJf,aii''^o) is shown 
in Fig. ^ (right). It reveals a smooth mass dependence of the scaling optimal overlap HF. 

Finally we remark that we can also insert a massive Dq , and make the resulting overlap 
fermion massive again. As an example, we consider some massive TP-HFs with varying 
mass parameters M, which are all converted to a massive overlap fermion by adding m = 1. 
Here the static masses mo come out differently depending on M. For a direct comparison 
we show in Fig. ^ (right) the difference between the energy E{pi) and the continuum 
energy Econt{pi)- Again the vicinity of M = 1 looks very good, but — as we mentioned 
before — that fermion might be too close to zero- or multi-species (doubling) to be useful 
in gauge theory. 
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